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Supernova remnants (SNRs) have been regarded for many decades as the sources of Galactic cosmic rays (CRs) 
up to a few PeV. However, only with the advent of Fermi-LAT it has been possible to detect — at least in some 
SNRs — 7-rays whose origin is unequivocally hadronic, namely due to the decay of neutral pions produced 
by collisions between relativistic nuclei and the background plasma. When coupled with observations in other 
bands (from radio to TeV 7-rays), Fermi-LAT data present evidence for CR spectra significantly steeper than the 
standard prediction of diffusive shock acceleration, forcing us to rethink our theoretical understanding of efficient 
particle energization at strong shocks. We outline how, by including the effects of CR-triggered magnetic field 
amplification, it is possible to reconcile non-linear models of diffusive shock acceleration with 7-ray observations, 
in particular providing a successful application of such a theory to Tycho's SNR. Finally, we show how kinetic 
simulations can investigate the microphysics of the non-linear coupling of accelerated particles and magnetic 
fields, probing from first principles the efficiency of the Fermi mechanism at strong shocks 



1. THE SNR PARADIGM 



2. EVIDENCE OF STEEP SPECTRA 



Even before the launch of the Fermi satellite, several 
supernova remnants (SNRs) have been detected in the 
7-rays by Cherenkov telescopes, attesting to the pres- 
ence of multi-TeV particles at their blast waves. This 
evidence supports the SNR paradigm for the origin of 
Galactic cosmic rays (CRs), which relies on the fact 
that: i) SNRs can accoun t for the proper energetics 
Baade and Zwicky||1934| ; ii) non-relativistic, strong 



shocks are expected to accelerate particles distributed 
in energy as an universal power-law , oc E~ 2 [e.g. |Bell| 
1978I iBlandford and Ostriker||1978|. In SNRs signa- 



tures of accelerated electrons have been found since 
the '50s, in the shape of radio synchrotron emission; 
however, the presence of relativistic nuclei can be di- 
rectly probed only in the 7-rays resulting from the 
decay of 7Tos produced in nuclear collisions between 
accelerated ions and the background plasma 

In the TeV band, the 7-ray spectrum almost in- 
variably shows a high-energy cut-off, which reflects 
the cut-off of the distribution of accelerated particles. 
With only this information available, it is impossible 
to disentangle whether the emission is hadronic, i.e., 
due to neutral pion decay, or leptonic, i.e., due to the 
inverse- Compt on scattering of relativistic electrons on 
some photon background. 

The biggest difference between the two scenarios 
above is the slope of the photon energy distribu- 
tion: while pion decay produces a spectrum par- 
allel to the proton one, inverse-Compton produces 
a significantly flatter spectrum. Since acceleration 
is rigidity-dependent, at energies where synchrotron 
losses are negligible, electron and ion spectra are par- 
allel (ex E~ q ). As a consequence, the expected pho- 
ton spectrum is oc ^ _ ( < ?+ 1 )/ 2 (ex v~ q ) in the leptonic 
(hadronic) scenario. The interesting range of photon 
energies that allows us to test the very nature of the 
7-ray origin falls exactly in the Fermi-LAT band. 



2.1. 7-ray emission from SNRs 

From a collection of the spectral indexes of SNRs 
detected in the GeV and/or in the TeV ba nd |Caprioli 
|2Q11| |Mandelartz and Becker Tjus||2Q13] , it is rather 
clear that the leptonic scenario is often disfavored: ex- 
cept RX Jl 713. 7-3946 (and, possibly, Vela Jr.), all of 
the sources show spectra steeper than E~ 2 . Quite in- 
terestingly, when a SNR is detected both in the GeV 
and in the TeV band, its spectrum can usually be fit- 
ted with a single power-law, suggesting a coherent ori- 
gin for the 7-ray emission. The alternative is to fine- 
tune the leptonic contributions of inverse-Compton 
and relativistic bremsstrahlung, as it has been pro- 
posed, e.g., for Cas A Abdo et aL|2010 



The photon spectral slope typically reads q 1 = 
2.2 — 2.5; such 7-ray spectra, to be explained as due to 
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Figure 1: Spect ral slopes, g 7 , o f 7-ray-bright SNRs. 
From table 1 in | Caprioli [2011 (see references therein), 
updated with data of W41, Puppis A and Cygnus Loop. 
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inverse- Compton, would require an electron distribu- 
tion as steep as E~ 3 - 5 — E~ 4 , at odds with radio/X-ray 
observations, and with any rigidity-dependent accel- 
eration mechanism at strong shocks. 

2.2. SNRs and molecular clouds 

Some 7-ray-bright SNRs, however, may belong to 
a distinct class of sources because of their interaction 
with mole cular clouds (MCs) [see, e.g., Castro and 
Slane|2010| . Some SNRs (W44, W49B, W28N and IC 
443, for instance), show a peculiar cut-off around 3- 
10 GeV, possibly related to the interaction of escaping 
CRs with MCs as well. The correlation of 7-rays with 
high gas densities is a strong hint that the emission 
is of hadronic nature. Very recently |Ackermann et al] 
|2Q13| confirmed this hypothesis by detecting the char- 
acte ristic pion-decay feat ure in IC 443 and W44 [see 
also |Giuliani et aL 2011 for AGILE measurements]. 
Nevertheless, one must bear in mind that 1) MC s may 
shine also because of diffuse Galactic CRs [e.g |Aha-| 
ronian et al.| 2006]; 2) modeling the expected spec- 



trum from a MC is not straightforward, because of the 
poorly-understood diffusion of CRs close to sources 



and in very dense environments iGabici et al. 2007 



SNRs close to, or interacting with, MCs show hadronic 
emission, but might not be probing freshly accelerated 
particles. Certainly, they are very hard to model in 
a self-consistent fashion, also because of the lack of 
emission in other bands of the spectrum. 

2.3. Anisotropy of diffuse Galactic CRs 

There is another observational evidence that the 
spectrum injected in the Galaxy by SNRs must be 
steeper than E~ 2 . The primary to secondary ratios 
and the anisotropy in the CR fluxes observed at Earth 
can be accounted for only if the residence time in the 
Milky Way has a mild dependence on the CR rigid- 
ity: oc (E/Z)~ s , with 6 ~ 1/ 3 [e.g 



2012, Blasi and Amato 2012 . Since S represents the 



Vladimirov et al. 



correction to the injection spectrum due to propaga- 
tion in the Galaxy, the constraint q + 5 ~ 2.75 implies 
that SNRs inject CRs with spectra steeper than E~ 2 . 
Differential escape of CRs can steepen the instanta- 
neous spectra probed in the 7-rays by 0.1 — 0.2 in slope 
Caprioli et al.pOlO 2011 , but the most straightfor- 



ward scenario points towards particles being acceler- 
ated with spectra oc E~ 23 — E~ 2A . 



3. THEORY VS OBSERVATIONS 

Steep spectra are also a challenge to the current the- 
ory of diffusive shock acceleration (DSA). The general 
consensus among non-linear models of DSA, namely 
models that account for the dynamical back-reaction 
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Figure 2: Top panel: downstream magnetic field, Bd, and 
CR acceleration efficiency, P cr , as a function of the 
fraction of injected particles, n. Bottom panel: slope q 
obtained with the scattering center velocity calculated in 
the amplified va(SB), or in the background field 
(Standard Non-Linear DSA). The values a re for a SNR in 
the early Sedov stage [about 3000yr, from Ca prioli|2012| . 



of CR pressure and energy, has traditionally been that 
the more efficient the particle acceleration, the flatter 
the CR spectrum [see |Malkov and O'C. Drury||2001[ 
for a review]. E~ 2 is the steepest spectrum expected 
in standard non-linear DSA, and it corresponds to a 
very inefficient particle acceleration (less than 1%), at 
odds with the requirement of the SNR paradigm (10- 
30%). This is not a small correction, but a systematic 
discrepancy: if acceleration is efficient, the canonical 
theory predicts the CR spectrum to be invariably flat- 
ter than in the test-particle limit. In the next section 
we outline a possible way of reconciling the theory of 
efficient DSA with the steep spectra observed in SNRs. 



3.1. Magnetic field amplification 

SNR blast waves are associated with enhancements 
of the magnetic field by factors of 10 — 100 with re- 
spect to the mean Galactic value, as it follows from 
several observational facts: 1) young SNRs show non- 
thermal X-ray rims, whose thickness can be explained 
as the synchrotron loss-length of multi-TeV electrons 
radia ting in magnetic fields as large as 100 — 500/iG 
[e.g.,|Volk et al.||2005[ [Parizot et al.||2006|; 2) X-ray 



variability on a time-scale of a few years is consis- 
tent with — localized — mG magnetic fields Uchiyama 
et al. 2007] ; 3) fitting radio-to-X-ray data typically 
shows electrons to be synchrotron-cooled above a crit- 
ical energy, at which the loss-time equals the age of 
the remnant. In Tycho, for instance, this constraint 
returns an average post-shock field of about 200/iG 



Morlino and Caprioli||2012|; 4) the lack of detection, 



within Chandra's resolution, of X-rays ahead of the 
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shock of SN1006 suggests that the field is amp lified in 
a CR- induced precursor |Morlino et al. 2010] . Mag- 
netic field amplification is also important to effectively 
scatter particles back and forth across the shock, po- 
tentially allowing SNRs to accelerate protons up to a 
few PeV, and heavier ions to larger energies, propor- 
tionally to their charge. Self-generated fields may even 
be dynamically relevant, in that the magnetic pressure 
may become comparable to, or even larger than, the 
plasma pressure, hence reducing the pre-sho ck mag- 
net osonic Mach number Caprioli et al.||2008 . 



The most intriguing aspect of magnetic field am- 
plification, though, is that these strong fields are ex- 
pected to be generated by plasma instabilities driven 



by th e same accelerated particles [see Bykov et al. 



2012, for a recent review]. If the growth rate is large 
enough to allow the turbulence to become nonlinear, 
as required to amplify the field, the local Alfven ve- 
locity may increase significantly. If we take the Alfven 
velocity va(SB) as a proxy for the velocity of the 
magnetic irregularities particles scatter against, the 
CR diffusion- convection equation has to be accord- 
ingly modified, eventually leading to a modification 
of the CR spectrum. The very reason behind this 
phenomenon is that the CR spectral slope is actually 
determined by the velocity jump r = v u /vd felt by 
particles diffusing between upstream and downstream 
(or viceversa). For strong gaseous shocks, r = 4 and 
Q — ^rrf = 2; if relativist ic CRs are present, r > 4 and 
q < 2. However, if va is not negligible with respect 
to v u , then we may have r = Vu ~^ A < 4 and q > 2, 
at least in the sim ple case of excitation of resonant 
Alfenic modes 



[e.g. 



Bell 1978 



where self-generated 
waves upstream propagate against the fluid (i.e., in 
the direction opposite to the CR gradient), and down- 
stream are almost isotropic (vA,d ~ 0). 

If acceleration is efficient, the CR pressure tends to 
make the spectrum flatter, while field amplification 
may provide an opposite backreaction that steepens 
the spectrum. In order to see which effect dominates 
over the other, a self-consistent model is required. Us- 
ing a semi-analytical approach to non-linear DSA, it is 
possible to show that the more efficient the CR accel- 
eration, the more effective the magnetic field amplifi- 
cation, and the steeper the spectra of the accelerated 
particles 



Caprioli 2012 



Such a trend is showed in 
Fig.[2j where CR pressure, downstream magnetic field, 
and spectral slope are plotted against the fraction of 
particles injected into the acceleration mechanism, 77, 
for a SNR in its early Sedov stage [see Caprioli||2012 
for more details]. For rj > 10 -4 the acceleration effi- 
ciency P cr (expressed as the pressure in CRs over the 
bulk pressure) saturates to 10 — 20%, but the spectrum 
becomes steeper and steeper. The non-linear feedback 
of the amplified fields may account, at the same time, 
for efficient acceleration, for the amplified magnetic 
fields (a few hundred /iG), and for the spectral slopes 
inferred from 7-ray observations (q = 2.1 — 2.5). 
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Figure 3: Tycho's broadband emission. Curves show the 
calculated synchrotron, inverse Compton, bremsstrahlung 
and p ion decay emission. Refer to |Morlino and Caprioli| 



2012 



for the experimental data and for more details. 



3.2. Tycho's SNR 

Tycho is probably the best SNR where to study 
DSA: it corresponds to the historical SN1572 (so that 
its age is known), and it is the remnant of a type-la 
explosion, which puts a constraint on the SN energy 
(« 10 51 erg) and on the ejecta mass (about one solar 
mass); both kinematic and echo-light measurements 
set its distance within 3 and 4 kpc, implying a cir- 
cumstellar density of ~ 0.3 protons cm -3 . Also, de- 
tailed maps in radio and X-rays are available, and the 
shell presents a rather clean spherical symmetry. All 
of this information significantly reduces the number of 
free parameters, making modeling Tycho's broadband 
emission an over-constrained problem. There are two 
parameters that need to be tuned: the fraction of pro- 
tons injected into the acceleration mechanism, and the 
electron-to-proton ratio at relativist ic energies. The 
maximum proton energy can be estimated by equat- 
ing the acceleration time with the age of the system, 
which is also consistent with requiring the diffusion 
length of the highest-energy particles to be about one 
tenth of the SNR radius. Given a reliable analytical 
description of the SNR evolution Truelove and Mc 
Kee|19 99 , a self-consistent treatment of the shock dy- 
namics including CR and magnetic contributions, and 
a model for magnetic field amplification, it is possible 
to work out the expected mult i- wavelength, radially- 
resolved, map of the non-thermal emission with two 
free parameters only |Morlino and Caprioli|2012| . 

The less solid parts of the calculation are the es- 
timate of the growth rate of the resonant streaming 
instability, extrapolated from the quasi-linear regime, 
and the assumption that the Alfven velocity traces 
the velocity of the scattering centers. Also, we assume 
Bohm diffusion, in the amplified field. Strictly speak- 
ing, the Bohm limit corresponds to resonant scat- 
tering onto self-generated Alfven waves of amplitude 
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SB/B ~ 1, and it is achieved when the CR energy 
distribution is cx E~ 2 . 



Recently, iBerezhko et al. 20131 criticized the adop 



tion of a Bohm-like diffusion coefficient when the CR 
spectrum is steeper than E~ 2 , arguing that in this 
case diffusion should be less effective, thereby reduc- 
ing the maximum energy achievable at any given time. 
Firstly, the proton spectrum we obtained retains a 
hint of the curvature typical of CR- modified shocks, 
flattening at the highest energies, as pointed out in 
fig. 4 of |Morlino and Caprioli| |2Q12| . Secondly, and 
most importantly, when SB/B ^> 1 Bohm diffusion is 
not more than a heuristic assumption: the normaliza- 
tion of the diffusion coefficient, adopted by Berezhko 
|et al.|[2Q13] as well, must be a factor of SB/B ~ 30— 50 
smaller than the one predicted by the quasi-linear the- 
ory to allow CR acceleration up to a few hundred 
TeV. In this respect, a deviation of less than 10% 
in the energy dependence of the diffusion coefficient 
seems indeed acceptable. The main criticism than 
can be raised against all of the models proposed is 
that the magnetic field must be amplified even far up- 
stream, on the diffusion scale of the highest-energy 
particles. Resonant streaming instability is unlikely 
to be relevant there, and magnetic field amplification 
has probably to rely on instabilities triggered by es- 
caping particles, like the non-resonant hybrid insta- 
bility Bell 2004 . If this is the case, however, it is 
still not clear whether the excited modes, which have 
short-wavelengths and are non-resonant with ions in 
terms of polarization, can scatter CRs effectively. 

A crucial difference between our model and the one 
proposed by 



Berezhko et al. 



2013 



is that in ours 
the magnetic field is calculated according to a sim- 
ple but physically- motivated recipe, and not put by 
hand. Therefore, the spectral steepening induced by 
the enhanced velocity of the scattering centers is not 
arbitrary, but it is consistently related to the resulting 
CR acceleration efficiency and to the inferred mag- 
netic fields. Conversely, IBerezhko et al. 20131 explain 



Tycho's steep 7-ray spectrum by introducing an addi- 
tional contribution of CR protons accelerated at sev- 
eral low-Mach-number reverse shocks produced by the 
interaction of the blast wave with small, dense clouds 
filling the circumstellar medium. This scenario, how- 
ever, requires quite a fine-tuning of the properties of 
these weak shocks in terms of densities, temperatures, 
velocities, magnetic fields, and electron synchrotron 
emissivity to match the 7-ray emission without violat- 
ing any other observable, which the authors fit with 
the " canonical" CR population only Volk et al.|2008 



It is quite unlikely for such a conspiracy to work also 
in many other SNRs, while steep spectra are ubiqui- 
tous (Fig. [I]). In this respect, it seems more natural 
to conceive the hypothesis that CRs may be acceler- 
ated with spectra steeper than those predicted by the 
standard non-linear DSA theory. 




Figure 4: Evolution of the post-shock energy spectrum in 
a hybrid simulation of a parallel M = 20 shock. Time is 
units of (ion cyclotron frequency) -1 , and energy is in 
E s h = ^Vshi with V s h the shock velocity. The spectral 
slope matches the DSA prediction for strong shocks. 



4. HYBRID SIMULATIONS OF 
COLLISIONLESS SHOCKS 

As it is clear also from the discussion above, any the- 
oretical model aiming to account for the many obser- 
vational constraints requires a detailed understanding 
of the plasma instabilities triggered by the accelerated 
particles, including the nature of the fastest-growing 
modes, their saturation, and their effects in scattering 
CRs. Also, particle injection has not been understood 
from first-principles, yet, and any non-linear approach 
to shock acceleration has to model this crucial ingre- 
dient through some phenomenological recipe. Finally, 
we need to prove that Fermi acceleration may be as 
efficient as inferred, and that it can sustain itself by 
generating enough magnetic turbulence to allow the 
maximum energy to increase with time. 

These crucial pieces of physics involve very small 
plasma scales, and must be tackled with numerical 
approaches. Particle-in-cell techniques can capture 
the intrinsically nonlinear nature of the interplay be- 
tween accelerated particles and the electromagnetic 
field: single-particle trajectories are calculated accord- 
ing to the Lorentz force produced by self-consistent 
currents. Nevertheless, PIC simulations are computa- 
tionally very expensive. To mitigate this limitation, 
the so-called hybrid limit can be used: (massless) elec- 
trons are taken as a neutralizing fluid, while ions are 

for a review!; 



treated kinetically [see Lipatov 



2002, 



this approach still captures both the shock dynam- 
ics and the ion anisotropy, which drives the plasma 
instabilities relevant for magnetic field amplification. 

Here we summarize some results recently obtained 
with the newtonian (i.e., non-relativistic) code dHy- 
brid |Gargate et al.| |2007 . We consider a parallel 



shock (background field Bq parallel to the shock veloc- 
ity V 3 h) with equal sonic and Alfvenic Mach numbers, 
M c± 20. The computational box measures 10 5 x 10 2 in 
units of ion skin depths c/uo pi with u p the ion plasma 
frequency, and 2 cells per ion skin depth are used; de- 
spite the box being 2D, three components of particle 
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CR-driven filamentation instability: the streaming of 
accelerated particles upstream of the shock produces 
cavities of rarefied plasma and low-magnetic field, sur- 
rounded byan^tofd^nse filaments with strong fields 
[see also Reville and Bell||2Q12 . The cavitation of the 
upstream, on scales comparable with the gyroradii of 
the accelerated ions, produces a significant increase 
of the magnetic field and of the local Alfven veloc- 
ity (see Fig. [5]); this effect may reveal itself crucial 
for enhancing particle diffusion. Moreover, the ad- 
vection of these cavities through the shock produces 
bubbles subject to (impulsive) Rayleigh-Taylor insta- 
bility, which triggers turbulent downstream motions 
enhancing the magnetic field amplification. Some 
knots in Fig. [5] carry a magnetic field ~ 40 times 
larger than £?o, and preliminary runs at M = 50 show 
that the maximum field can easily become as large as 
^ 100£>o. This localized strong fields may be consis- 
tent, for instance, with the X-ray variability reported 
by Uchiyama et al. 2007 . The downstream bubbles 



are filled by plumes of hot plasma, which stretch the 
magnetic field along the shock direction of propaga- 
tion (Fig. [5j density, B x , and B tot plots). These co- 
herent, elongated structures, which are predicted to 
develop only where the shock is parallel, may have 
an observational counterpart in the pattern of radial 
X-ray stripes detected in Tycho Eriksen et al.poTT 
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Figure 5: Density, parallel (B x ), transverse (B y ,B z ) and 
total (Btot) magnetic field, and Alfven velocity va for a 
2D simulation of a M = 30 parallel shock, at t — SOOco^T 1 . 
Quantities are normalized to their unperturbed values. A 
filament with B to t ~ 20^ and knots with B tot 40B 
are clearly visible [Caprioli and Spitkovsky|2013| . 



velocities and electromagnetic fields are accounted for. 

Fig. [4] shows the evolution of the post-shock ion 
energy distribution: the spectrum develops a non- 
thermal power-law tail, whose extension increases 
with time. Very interestingly, the spectral slope of 
the accelerated particles is the one predicted by Fermi 
acceleration at strong shocks, namely f(p) oc p~ 4 in 
momentum space, which translates into E~ 1,b for non- 
relativistic particles (p 2 = 2mE). If the spectrum ex- 
tended into relativistic energies, p oc E, and the stan- 
dard f(E) oc E~ 2 would be recovered. The total en- 
ergy in the non-thermal tail is about 10% of the shock 
bulk energy, attesting to the ability of parallel shocks 
to efficiently accelerate ions [Caprioli and Spitkovsky, 
in preparation]. Kinetic simulations do not allow to 
test all of the nonlinear effects mentioned above, yet, 
but they allow us to self-consistently study the devel- 
opment of the instabilities induced by energetic par- 
ticles. U nprecedentedly- large 2D a nd 3D hybrid sim- 
ulations |Caprioli and Spitkovsky 2013 



5. CONCLUSIONS 

More than a century after the discovery of CRs, and 
many decades after the discovery that electrons are 
accelerated in SNRs, we are finally getting direct evi- 
dence of the acceleration of ions in these environments. 
GeV observations by Fermi-LAT are complementing 
TeV data obtained with Cherenkov telescopes, open- 
ing an unprecedented window on the non-thermal 
emission of SNRs. In some cases th e 7-rays observed 
are unequivocally of hadronic origin Ackermann et al. 



2013 , but extracting quantitative constraints on the 



mechanisms responsible for ion acceleration requires a 
non-trivial theoretical modeling. When the emission 
comes from MCs illuminated by nearby SNRs, there is 
a degeneracy between the injection spectrum and the 
correction due to the poorly-known diffusion around 
the sources. Also, when shocks are directly interact- 
ing with dense, partially- neutral gas, the dynamics 
is strongly pertur bed by charge-exchange processes 
Blasi et aT| 2012 . Modeling these shocks requires 
a kinetic description of the neutral distribution, but 
the effort is rewarded with an additional diagnostics, 
based on the study of Ha emission , of ion tempera- 
ture and CR acceleration efficiency Ghavamian et al. 



20071 |Raymond efaLl[20TT| |Morlino et aTf2 012 . 



showed that 

high-Mach number parallel shocks are prone to the 



The most straightforward way to probe CR accel- 
eration at SNRs is to look at SNRs non interacting 
with MCs [which include Tycho: see Tian and Leahy 
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2Q11| , whose steep 7-ray spectrum challenges the cur- 
rent models of efficient particle acceleration at shocks. 
A crucial open problem is to understand in greater 
detail the mechanisms that amplify the magnetic field 
up to the inferred values, mechanisms that may be 
important also for regulating the spectral slopes of ac- 
celerated particles. Kinetic simulations play a pivotal 
role in investigating the complex particle-field inter- 
play from first-principles. At the core of the accel- 
eration process, however, always lies the pioneering 
idea envisioned by Enrico Fermi, and namely that a 
particle can increase its ener gy by b eing repeatedly re- 
flected in head-on collisions |Fermi|p~949 . More than 
60 years later, a satellite named after him is providing 
us with very strong evidence that SNRs are indeed the 
sources of Galactic CRs. 



Acknowledgments 

The author wishes to thank A. Spitkovsky and 
G. Morlino for their collaboration, and L. Gar- 
gate for his help with dHybrid. This research was 
supported by NSF grant AST-0807381 and NASA 
grants NNX09AT95G and NNX10A039G. Simulations 
were performed on the computational resources sup- 
ported by the PICSciE-OIT TIGRESS High Perfor- 
mance Computing Center and Visualization Labora- 
tory. This research also used the resources of the Na- 
tional Energy Research Scientific Computing Center, 
which is supported by the Office of Science of the U.S. 
Department of Energy under Contract No. DE-AC02- 
05CH11231, and Teragrid/XSEDE's Ranger under 
contract No. TG-AST100035. 



References 

W. Baade and F. Zwicky, Proceedings of the National 

Academy of Science 20, 259 (1934). 
A. R. Bell, MNRAS 182, 147 (1978). 
R. D. Blandford and J. P. Ostriker, ApJL 221, L29 

(1978). 

D. Caprioli, JCAP 5, 26 (2011), 1103.2624. 

M. Mandelartz and J. Becker Tjus, ArXiv e-prints 

(2013), 1301.2437. 
A. A. Abdo et al., ApJL 710, L92 (2010), 1001.1419. 
D. Castro and P. Slane, ApJ 717, 372 (2010), 

1002.2738. 

M. Ackermann et al., Science 339, 807 (2013). 
A. Giuliani et al., ApJL 742, L30 (2011), 1111.4868. 
F. Aharonian et al., Nature 439, 695 (2006), astro- 
ph/0603021. 

S. Gabici, F. A. Aharonian, and P. Blasi, APSS 309, 
365 (2007), astro-ph/0610032. 



A. E. Vladimirov, G. Johannesson, I. V. Moskalenko, 
and T. A. Porter, ApJ 752, 68 (2012), 1108.1023. 

P. Blasi and E. Amato, JCAP 1, 011 (2012), 
1105.4529. 

D. Caprioli, E. Amato, and P. Blasi, APh 33, 160 

(2010) , 0912.2964. 

D. Caprioli, P. Blasi, and E. Amato, APh 34, 447 

(2011) , 1007.1925. 

M. A. Malkov and L. O'C. Drury, Rep. Prog. Phys. 
64, 429 (2001). 

D. Caprioli, JCAP 7, 038 (2012), 1206.1360. 

H. J. Volk, E. G. Berezhko, and L. T. Ksenofontov, 
A&A 433, 229 (2005), astro-ph/0409453. 

E. Parizot, A. Marcowith, J. Ballet, and Y. A. Gal- 
lant, A&A 453, 387 (2006), astro-ph/0603723. 

Y. Uchiyama, F. A. Aharonian, T. Tanaka, T. Taka- 
hashi, and Y. Maeda, Nature 449, 576 (2007). 

G. Morlino and D. Caprioli, A&A 538, A81 (2012), 
1105.6342. 

G. Morlino, E. Amato, P. Blasi, and D. Caprioli, MN- 
RAS 405, L21 (2010), 0912.2972. 

D. Caprioli, P. Blasi, E. Amato, and M. Vietri, ApJL 
679, L139 (2008), 0804.2884. 

A. M. Bykov, D. C. Ellison, and M. Renaud, Space 

Sci. Rev. 166, 71 (2012), 1105.0130. 
J. K. Truelove and C. F. Mc Kee, ApJ Supplement 

Series 120, 299 (1999). 

E. G. Berezhko, L. T. Ksenofontov, and H. J. Volk, 
ApJ 763, 14 (2013), 1211.5398. 

A. R. Bell, MNRAS 353, 550 (2004). 

H. J. Volk, E. G. Berezhko, and L. T. Ksenofontov, 
A&A 483, 529 (2008), 0803.1403. 

A. S. Lipatov, The hybrid multiscale simulation tech- 
nology: an introduction with application to astro- 
physical and laboratory plasmas (Berlin; New York: 
Springer, Scientific computation, 2002). 

L. Gargate, R. Bingham, R. A. Fonseca, and L. O. 
Silva, Comp. Phys. Commun. 176, 419 (2007), 
arXiv:physics /061 1 174. 

D. Caprioli and A. Spitkovsky, ApJL 765, L20 (2013), 
1211.6765. 

B. Reville and A. R. Bell, MNRAS 419, 2433 (2012), 
1109.5690. 

K. A. Eriksen et al., ApJL 728, L28+ (2011), 
1101.1454. 

P. Blasi, G. Morlino, R. Bandiera, E. Amato, and 

D. Caprioli, ApJ 755, 121 (2012), 1202.3080. 
P. Ghavamian, J. M. Laming, and C. E. Rakowski, 

ApJL 654, L69 (2007), astro-ph/0611306. 
J. C. Raymond, J. Vink, E. A. Helder, and A. de Laat, 

ApJL 731, L14 (2011), 1103.3211. 
G. Morlino, P. Blasi, R. Bandiera, E. Amato, and 

D. Caprioli, ArXiv e-prints (2012), 1211.6148. 
W. W. Tian and D. A. Leahy, ApJL 729, L15 (2011), 

1012.5673. 

E. Fermi, Physical Review 75, 1169 (1949). 



eConf C121028 



